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We present experimental and modeling studies of UV nanosecond pulsed laser desorption and 
ablation of (111) bulk silicon. The results involve a new approach to the analysis of plume formation 
dynamics under high-energy photon irradiation of the semiconductor surface. Non-thermal, photo- 
induced desorption has been observed at low laser fluence, well below the melting threshold. Under 
ablation conditions, the non-thermal ions have also a high concentration. The origin of these ions 
is discussed on the basis of electronic excitation of Si surface states associated with the Coulomb 
explosion mechanism. We present a model describing dynamics of silicon target excitation, heating 
and charge-carrier transport. 
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I. INTRODUCTION 



(the Coulomb explosion mechanism). 
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Laser ablation technique has been widely used for thin 
film deposition and nanoparticles synthesis. For these ap- 
plications, understanding of dynamics and mechanisms 
of particle desorption and ablation from the irradiated 
surfaces is of crucial importance. For nanosecond laser 
pulses, it is generally accepted that the ablation mech- 
anism of normal vaporization gives way to phase explo- 
sion with increasing laser fluence Q. At low laser flu- 
ences near ablation threshold, electronic mechanism of 
high-energy ion emission has been proven to play a role 
in initiating the ablation process, mainly for dielectrics 
and semiconductors 0, ■ For ultrashort laser pulses, 
the electronic mechanisms of desorption Q and ablation 

|ij for both dielectric and semiconductor targets are 
studied extensively, whereas for nanosecond pulses there 
is a lack of both the experimental and theoretical analy- 
ses. 

In this paper, we present the results of the experi- 
mental and theoretical studies of ion ejection from the 
Si targets induced by UV nanosecond laser pulses in a 
wide range of laser fluences in the regimes from those 
well below melting threshold up to developed ablation. 
Non-thermal ions have been detected at very low laser 
fluence. The theoretical study of ion ejection involves 
an analytical analysis and numerical modeling of charge 
transport in the laser-irradiated target. We attribute the 
high-energy ion emission to the generation of a strong 
electric field due to the electron photoemission process 
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II. EXPERIMENTAL 

Clean (111) Si surface was passivated by hydrogen in a 
usual chemical way Q ■ Ablation was carried out in a high 
vacuum chamber (pressure < 1CU 8 Torr) at an incidence 
angle of 45 ° using an ArF excimer laser (hv = 6.4 eV, 
15 ns pulse duration at FWHM). The set of masks was 
used to select the homogeneous part of the laser beam. 
The laser spot size on the target and laser fluence (Fq) 
were 0.5 mm 2 and 0.01 — 1 J/cm 2 , respectively. The 
target was rotated/translated during measurements to 
avoid cratering. The expansion dynamics and origin of 
the desorbed particles were analyzed by reflectron time- 
of-flight (TOF) mass spectrometer. 

The most abundant desorbed species, Si + , are ob- 
served already at Fo ~ 0.2 J/cm 2 , well below melt- 
ing threshold determined by time-resolved reflectivity 
measurements to be ~ 0.4 J/cm 2 (see Q and refer- 
ences therein) that is consistent with the other studies 
0, [^} . Neutral monatomic Si particles are detected 
only at Fo > 0.8 J/cm 2 . Attempts to detect charged 
or neutral clusters, excepting silicon ions dimmer [llj, 
have been unsuccessful. Figure [I] shows the typical TOF 
spectra of Si + at different fluences. At low laser fluence 
(0.26 J/cm 2 on the Fig. ^ the TOF spectrum consists 
of an only population of Si + ions. The spatial distri- 
bution of this kind of ions is narrow, strongly peaked 
relatively to the target normal and, up to the melting 
threshold, the kinetic energy (~ 4.8 — 5 eV) is a weak 
function of laser fluence that is typical for non-thermal 
desorption / ablation 0, 0] • 

Increasing laser fluence above the melting threshold 
leads to broadening of Si + TOF distribution (see Fig. ^ 
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FIG. 1: Time-of- flight spectra of Si + ions for different laser 
fluences: 

open circles - 0.26 J/cm 2 intensity (multiplied by 100); 
dark triangles - 0.42 J/cm 2 ; 
open triangles - 1.02 J/cm 2 . 

The target - mass spectrometer distance was 126 mm. 
Inset shows charge of silicon target versus laser fluence. 
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FIG. 2: Experimental data on Si + yield vs. laser fluence. 
0.2 J/cm 2 corresponds to the ion desorption threshold. Up 
to ~0.4 J/cm 2 that is around melting threshold, ion yield be- 
haves as Fq 0A whereas at higher fluence dependence changes 
to F$- 97 . Solid curve corresponds to Eq. © with Foth = 0.2 
J/cm 2 . 



fluence 0.42 J/cm 2 ) and to formation (from 0.7 — 0.8 
J/cm 2 ) of a well-pronounced low-energy distribution [l3j . 
The appearance of the low-energy Si + concurs with the 
development of Si surface melting. The intensity of this 
population rapidly increases with fluence and becomes 
dominant at Fq > 1 J/cm 2 . The abundance of Si + 
species versus laser fluence (Fig. [21 experimental points) 
shows an abrupt change in the behavior of the exper- 
imental data that can be treated as the change of the 
ejection regime. The first regime, corresponding to the 
generation of the fast population, exhibits a very strong 
non-linear dependence on fluence, -Fq 10 ' 4 , that is a typi- 
cal feature of the multi-particle ejection process The 
second regime, starting after melting, shows near-linear 
Si + abundance variation with fluence. 

The expulsion of the non-thermal Si + ions is attributed 
to accumulation of a positive charge on the silicon sur- 
face as a result of electron photocmission. This charge 
has been evaluated by measuring a compensating current 
variation between the Si target and ground [l4| during 
laser pulse. The positive charge is accumulated already 
at ~ 0.001 J/cm 2 (see inset on Fig. [IJ, increasing near 
linearly up to a saturation level of (1 — 2) • 10 13 elementary 
charges per pulse at ~ 0.15 — 0.2 J/cm 2 . Exactly at these 
fluences we detect the first non-thermal ions. Both the 
target charge saturation and ions ejection clearly indicate 
that Si surface reaches the critical, threshold, conditions 
corresponding to dynamical equilibrium between primary 
particles ejection (photoelectrons) and desorption of the 
secondary species (positively charged ions). 

III. MODELING 

An assumption for the electron photoemission under 
the UV laser irradiation, when photon energy exceeds 
work function, reads as 

1 I(x,t) ( x \ 

where I(x,t) is the laser power of a Gaussian temporal 
shape: 

I(x,t) = (1 — R) ^yi^cxp \-4\n2 e- Q — , 

(2) 

t is laser pulse duration (FWHM), a a b and R are the 
absorption and reflection coefficients, Ipe is the electron 
escape depth, x is the distance from the target surface to- 
ward the bulk depth, hv is the energy of laser light quan- 
tum. The expression Q implies that the laser-generated 
electrons whose momentum component normal to and in 
the direction of the surface are immediately photoemit- 
ted from the surface and below surface region with an 
exponential decreasing within the bulk [f| . 
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An estimate of the number of electrons photoemitted through a unit surface area (Npp) is obtained by integrating 
the photoemission term over time and space that gives: 

a ab (l - R)F 
2hv{a ab + l pE ) 

Under the assumption that the target is unearthed and there is no electron supply from the radiation-free sides of 
the target, the electric field generated on the target surface can be estimated by using the Gauss law (the net positive 
charge of the target is equal to the charge of the photoemitted electrons) : 

eNpE ea ab {l - R)F 



E\ = — — / a (x)dx = = °° v — '— - f4") 

lx =° 2ee J KJ ' 2ee 4ee his(a ab + l PE ) ' K) 



where L is the target width and e is dielectric permittivity of bulk silicon. 

The threshold electric field necessary to be exceeded in order to break the atomic bonds in crystalline silicon can 
be estimated through the energy density of the electric field, w = eeoE 2 /2. The value W at = eeqE 2 jln falls at an 
atom in the crystal, where rt _1 is a volume occupied by atom. The binding energy of an atom on the target surface 
estimated from the latent heat of sublimation, A s „6 = 16115 J/g Q, is ~ 4.67 eV. Thus, the threshold electric field 
is of order E t h\x=o = \/2h. at n/ 'eeo or E t h ~ 2.65T0 10 V/m for crystalline silicon. Laser fluence above which the 
non-thermal ions can be observed is: 

a a b + ipg) , . 

roth = 7Z 7T\ ■ (Oj 

ea ab \l - R) 

With Ipe ~ 10 A it gives ~ 0.13 J/cm 2 that is in good agreement with the experimental observations. The excess 
of positive charge in the target for generation of the threshold electric field is eN Ith = 2 ee SpE th \ x=0 or ~ 1.8 • 10 13 
electrons escaped from the irradiation spot size Sp of 0.5 mm 2 that is in excellent agreement with the measured 
electric residue. 

A strongly charged target tends to take off the electrostatic stress, decreasing the electric field below its threshold 
value through ion expulsion. The number of ions thrown out from the target by the electrostatic force can be evaluated 
as 

AWk = N, - N Ith = QQb(1 9 :f (F °~y* (6) 

2hv{a ab + l PE ) 

or 1.4-10 14 (i 7 b — Foth) for our irradiation conditions (laser fluence is taken in J/cm 2 ). This dependence follows the 
general tendency observed experimentally (Fig. [21 solid line). Thus, a seeming change in the ablation regime from 
the strong non- linear to near- linear dependence of ion yield vs. laser fluence can be roughly described by a shifted 
linear dependence [Eq. |JBJ]. 

More sophisticated description can be reached by modeling the electron dynamics in the laser-irradiated semicon- 
ductor targets. Here we shall consider only the target charging effect without concerning the ion ejection process, both 
electrostatic and thermal. Our model is based on the continuity equations for electron and hole generation, including 
one-photon ionization, Auger recombination and the photoemission process according to Eq. JQ), and we incorporate 
the electric current J y (y = e, h) = |e| n y [iE — eDS7n y written in the drift-diffusion approach into the equations: 

dn e e dn e d dn e I(x, t) 2 

— — - H e n e [n h - n e ) - HeE— — D e — — = a ab — 7«e n h _ P E, (7) 

at ee a ox ox Ox nw 

dn h e dn h d dn h I(x,t) 2 

-57- + fJ-hrih (n h - n e ) + fihE— -z-D h —— = a ab — in e n h . (8) 

at eeo ox ox ox noj 

Here n, fi, D, 7 are the density, mobility and the diffusion and Auger recombination coefficients, respectively; indices 
e and h refer to the electrons and holes. The diffusion coefficients are calculated as D e = kBT e [i e /e and Dh = kpTe^hl e 
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with the carrier temperature T e . The photoemission and diffusion lead to charge separation in the target that results 
in the electric field generation described by the Poisson equation with the boundary condition according to the Gauss 
law. The energy equations for the electron and lattice subsystems degenerate into an only equation at fairly long 
laser pulses of ns time scale 0, 



c P P 



^L- u (t) — \ - —X— 
dt dx ) dx dx 



ow>(l - R)h{t) cxp(-a ab x) 



(9) 



The vaporization rate u(t) is defined under the assump- 
tion that the flow of vaporized material from the surface 
follows the Hcrtz-Knudscn equation and the vapor pres- 
sure above the vaporized surface can be estimated with 
the Clausius-Clapeyron equation: 



W 1 U p \2nkT s 



1/2 



exp 



f 

n 



i 



(10) 

Thermodynamic and optical properties of silicon (spe- 
cific heat c p , thermal conductivity A, latent heat of vapor- 
ization L, boiling temperature Tf,, absorption coefficient 
a a b, and reflection coefficient R) were taken from Ref. 
0- The initial and boundary conditions are 
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T(0,x)=T o , T(t,0) 



dT 

T ^ ^ ~d~x 



Lu(t) (11) 



with To being the initial temperature uniform across the 
target. The target was divided into an irregular grid, 
dense in the absorption region (cells of 5 A) and rarefying 
toward the target depth. We used an explicit scheme for 
solving the described system of equations. The time step 
was selected empirically to satisfy the numerical scheme 
stability and the approximation of the original equations. 



IV. RESULTS AND DISCUSSION 

The modeling results are presented in Figs. El El To 
the end of the laser pulse, the surface charge reaches al- 
most 60%, an unexpected value for silicon irradiated by 
low laser fluence. The densities of the electrons and holes 
are presented in Fig. |31 Considerable charging takes 
place in a surface zone of order of 2 nm wide (~ 5 — 6 
atomic monolayers). It is generally agreed that, because 
of efficient Auger recombination, a silicon target can not 
be ionized to a high ionization degree at moderate laser 
fluences. That is a major argument against the Coulomb 
explosion mechanism. Our modeling shows that a thin 
surface layer of the target reaches a high ionization de- 
gree, but therewith this layer occurs to be strongly de- 
pleted of electrons. Extremely high density of positive 
charge is gained in a few monolayers, implying the strong 
repulsion force between the ions. First of all, adatoms, 
which have lower binding with surface, have to be ejected 
from the surface, giving way to the generation of new 
adatoms. The process of ion expulsion will take place till 



FIG. 3: Calculated spatial profiles of the electron (dashed 
line) and ion densities (solid line) in the Si target irradiated 
with 0.2 J/cm 2 . 



all the layer of high positive charge is exploded or till the 
macroscopic electric field is reduced below the threshold 
value. Figure 01 illustrates that, at 0.2 J/cm 2 , the macro- 
scopic electric field exceeds the threshold value, as it was 
predicted analytically. The values of charging and the 
electric field in two external monolayers are so high, that 
their electrostatic explosion is inevitable. 

Let us consider the origin of high surface charging of 
a silicon target in detail. What is happening when a 
target subjected to laser irradiation loses the electrons 
due to photoemission? The dielectric breakdown takes 
place within a skin layer with an exponential decay to- 
ward the target depth. Because of photoemission, the 
target quasi-neutrality is broken. According to the Gauss 
law, the target, which starts to behave as a metal, aspires 
to accumulate the excess positive charge in its surface. 
The electric field in this thin surface layer is "negative" 
(directed to vacuum) that means that the non-emitted 
electrons are dragged toward the target depth to a re- 
gion with a lower field. Thus, the following scenario of 
the electrostatic mechanism for both ion desorption and 
ablation from the semiconductor targets (that may be a 
case also for dielectrics) can be proposed step by step: 
(i) target material breakdown and photoemission under 
laser irradiation; (ii) generation of the electric field which 
drags the electrons away from the surface layer; (iii) con- 
tinuing surface ionization together with suppression of 
Auger recombination; (iv) as a result, ionization degree 
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FIG. 4: Spatial distribution of the electric field in the Si 
target (F = 0.2 J/cm 2 ). 

up to 100% in a thin surface layer; (v) electrostatic disin- 
tegration of this layer. The most important consequence 



is that, in the charged outer layer, one-photon ionization 
is continuing during the laser pulse, whereas the Auger 
recombination is strongly suppressed because of lack of 
electrons. 



V. CONCLUSIONS 

Mechanisms of desorption and ablation of silicon un- 
der irradiation by the nanosecond laser pulses in a wide 
range of laser fluences have been studied both experimen- 
tally and theoretically. At low laser fluence the desorp- 
tion flux is mostly composed of high-energy ions with a 
narrow energy distribution that points to a non-thermal 
mechanism of their generation. Theoretical analysis has 
attributed the high-energetic ions to a strong electric 
field generated in an exterior target layer due to pho- 
tocmission that causes Coulomb explosion. More sophis- 
ticated two-tempcraturc modeling, which takes into ac- 
count electronic ion emission and the aspects of thermal 
vaporization, is under the progress. 
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